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Abstract 

We discuss how a L e — — L T flavour symmetry that is softly broken leads to keV 
sterile neutrinos, which are a prime candidate for Warm Dark Matter. This is to our 
knowledge the first model where flavour symmetries are applied simultaneously to 
active and sterile neutrinos explaining at the same time active neutrino properties 
and this peculiar Dark Matter scenario. The essential point is that different scales of 
the symmetry breaking and the symmetry preserving entries in the mass matrix lead 
to one right-handed neutrino which is nearly massless compared to the other two. 
Furthermore, we naturally predict vanishing #13 and maximal 623, while the correct 
value of #12 must come from the mixing of the charged leptons. We can furthermore 
predict an exact mass spectrum for the light neutrinos, which will be testable in the 
very near future. 
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1 Introduction 



The nature of Dark Matter (DM) is a longstanding problem of particle physics and astro- 
physics, see e.g. [1, 2, 3]. New elementary particles with suitable properties are predicted 
by various extensions of the SM and some of them are excellent candidates to make up the 
DM in the Universe. Any candidate does not only have to yield the right DM abundance, 
^dm^ 2 — 0.1 [ i], but it also has to be consistent with bounds coming from structure forma- 
tion provided by the Sloan Digital Sky Survey (SDSS) [5] and by the 2dF Galaxy Redshift 
Survey (2dF-GRS) [6]; see also the analysis of the so-called Lyman-a Forest [7]. These 
bounds rule out hot (relativistic) DM as the dominant component, which is why ordinary 
(active) neutrinos are excluded as Dark Matter candidate [4]. In addition active neutrinos 
are too light to explain the required DM abundance. A good alternative is that DM is cold 
(non- relativistic), which fits very well to the cosmological standard model (ACDM), which 
is in perfect agreement to observations [4]. 

A less explored possibility is so-called 'Warm Dark Matter', for which a single light sterile 
neutrino would be a natural candidate [8]. The advantage of such a scenario would be that 
right-handed neutrinos most likely exist anyhow in order to accommodate neutrino masses 
and one needs only some extra feature which explains why one sterile state is light in the 
multi keV range. Model-independent analyses of the data even seem to point into the 
direction of keV-masses for DM [9]. Such a scenario with three active plus one sterile light 
neutrino states fits also very well to the number of light neutrino species obtained from 
recent analyses of Big Bang Nucleosynthesis data [10]. A simple framework is to assume 
a certain parameter choice for the Standard Model (SM) extended by three right-handed 
neutrinos, which is called 'i^SM' [11]. This setting has been shown to be in agreement with 
cosmology [12, 13, 14, 15, 16, 17] and could also lead to detectable signals in laboratory 
experiments [18, 19, 20]. 

Another setting, which depends less on physics at even higher scale has recently been 
discussed by Bezrukov, Hettmansperger, and Lindner [21]. They studied right-handed 
neutrinos which couple to gauge extensions which naturally leads to the right DM abun- 
dance. In this context, it has been shown that keV sterile neutrino DM can be brought 
into agreement with current bounds in a type II seesaw framework, while type I scenarios 
are practically excluded. Apart from invoking a low scale seesaw mechanism, however, this 
paper did not provide an explanation why one sterile neutrino should have a keV-ish mass, 
while all the others are very heavy. One attempt to explain such a pattern was recently 
put forward by Kusenko, Takahashi, and Yanagida [22] in a scenario with extra dimensions 
leading to a 'split see-saw'. Here, we attempt to explain the required mass pattern by us- 
ing L e — — L T flavour symmetry in order to obtain keV sterile neutrino Dark Matter. 
Note that this would imply that the same flavour symmetry, which explain the masses and 
mixings of active neutrinos, would also affect the pattern of sterile neutrinos and explain 
in this way a natural DM scenario. Our study can be applied to the model of Ref. [21], 
but also to any other model that present a keV sterile neutrino as DM candidate. 
Our key idea is the following: Using L e — — L T , we can already explain hierarchical 
neutrino spectra with one neutrino out of three being massless. Soft breaking of this 
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symmetry brings in a completely new scale, much lower than the original one, which will 
determine the size of the induced mass of the neutrino that was massless before. Such a 
setup could very well point to scenarios with one sterile neutrino with a keV-mass, while 
the other two are heavier, which could be promising for cosmology as well as for laboratory 
detections [19]. The key point is that the same flavour symmetry can be responsible for 
similar structures in the heavy and in the light neutrino sector. In fact, all the cases studied 
by us predict a strong (inverted) hierarchy for the light neutrinos, which can be tested in 
the near future. 

Furthermore, the symmetry predicts bimaximal mixing coming from the neutrino sector. 
This involves a maximal mixing angle 6 , 12 , which is ruled out experimentally at more than 
6a. However, since there is also mixing coming from the charged leptons, it is possible to 
alleviate or even solve this problem. 

Our approach will be pragmatical in the sense that we simply consider an example of mixing 
matrix in the charged lepton sector that can lead to allowed values of neutrino mixing 
angles. We furthermore argue that it is possible to find charged lepton mass matrices 
which yield this mixing without contradicting L e — L yi — L T , in the sense that the largest 
elements of the mass matrix preserve the symmetry. Of course, one could try to go into 
some more detail here by having a closer look at the flavon potential and its breaking 
chain, in order to enforce this particular structure for the charged lepton mass matrix. 
This would, however, lead beyond the scope of this paper and would rather distract the 
reader from the simple statement we want to make, so we leave this point to further studies 
and focus on the neutrino part instead. 

The paper is organized as follows: In Sec. 2 we generally describe what mixing we need 
in the charged lepton sector in order to translate a bimaximal neutrino mixing into a 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix that is compatible with the experi- 
mental values. We then discuss our general framework and three different models in Sec. 3, 
before finally concluding in Sec. 4. 



2 The PMNS matrix 

Before we start with a discussion of the models under consideration, we shortly comment 
on the mixing that is predicted, and on how to bring it to agreement with experimental 
data. 

The PMNS matrix U is often parametrized as: 

^ c 12Cl3 S12C13 S13 \ 

U = -s 12 c 2 3 - C12S23S13 e lScp C12C23 - s 12 s 23 s 13 e iScp s 23 c 13 e iScp ■ diag(l, e ia , e i/3 ), 
\ S12S23 - C12C23S13 e l5cp -C12S23 - S12C23S13 e lScp C23C13 e tScp J 

(1) 

where Cij = cos^- and Sij = sin^-. The parameter 5qp is the Dirac CP- violating phase, 
while a and (3 are the two Majorana phases, absent in the case of Dirac neutrinos. The 
PMNS matrix is given by U = U\lA v , where U u is the matrix that diagonalizes the neutrino 
mass matrix, while Ul contributes to the diagonalization of the charged lepton mass matrix. 
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The PMNS mixing matrix should provide values of the neutrino mixing angles compatible 
with the 3a intervals allowed by the solar, atmospheric, reactor (KamLAND and CHOOZ), 
and accelerator (K2K and MINOS) experiments [23]: 

0.25 < sin 2 12 <0.37, 

0.36 < sin 2 # 23 <0.67, (2) 
sin 2 013 < 0.056. 



The allowed best-fit values and la errors of the mass square differences are given by 
Am^ = 7.65±g;|g x 10~ 5 eV 2 and lAmfJ = 2.401J];} 2 x 10~ 3 eV 2 . The sign of Am 2 31 is still 
unknown. It can be positive for normal mass ordering (mi < m 2 < m 3 ) or negative for 
inverted mass ordering (ma < m\ < m-z). No experimental information on the value of <5cp 
is present at the moment. 

A form of the mixing matrix which is frequently used and compatible with all current data 
is the so-called tri-bimaximal (TBM) mixing: 
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Details on the requirements for mass matrices in order to yield TBM can be found in 
Ref. [24]. Exact L e — — L T symmetry generically predicts a bimaximal form for the 



neutrino mixing matrix [25] 
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The PMNS matrix can deviate from the bimaximal structure reported above, depending 
on the actual form oilAi. We can parametrize the matrix Ul as a function of three mixing 
angles, in analogy to Eq. (1). For definiteness, we will focus on CP conservation and we will 
define Ajj = sin 6^-. Considering a hierarchical relation between A^-, of the type A12 = A, 
A13 ~ A 3 , A23 — A 2 , with A ~ 0.20 being the parameter that describes the deviation of 9 
from 7r/4 [26], we find 
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(5) 



It was shown in Ref. [26] that this form of Ul, combined with the bimaximal matrix of 
Eq. (4), could lead to Wpmns compatible with the experimental values: 

tan 2 12 ~ 1 - 2y/2X + 4A 4 - 2^2X 3 -> Q 
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The matrix IAl is associated to the diagonalization of the matrix AiiAi], with M.\ being 
the charged leptons mass matrix. The expression of Ul reported in Eq. (5) will diagonalize 
the matrix 

(m 2 . + m^X 2 m^X \ 
mlX ml , (9) 
m 2 T ) 

where we have neglected terms of order A 3 . In the next section we will show under which 
assumptions the charged lepton mass matrix M.iM.\ can be in agreement with (softly 
broken) L e — — L T symmetry. 

3 The L e — — L T flavour symmetry 

Let us define the symmetry T = L e — — L T . We extend the model [27] by a Higgs 
triplet A, in order to accommodate a type II seesaw, which is required in the context of 
the model presented in Ref. [21]. The particle content and the charge assignments of our 
model are given by: 
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where L aL = [y aL ,aj)} T , and a = e,/i,r. The scalars are defined as: 

A+/V2 A++ \ 
A -A+/V2J ' 

The right-handed neutrino fields are gauge singlets and we can thus write down the fol- 
lowing symmetry conserving Majorana mass term for the right-handed neutrinos: 

Cass = ~M l R 2 {NmfN m - M R 3 jN^fNtR + h.c. (10) 

The Dirac mass term, which links the left- and right-handed fields, is given by: 

Anass = -Y£lZL~0N 1R -Y£ 2 1^0N m -Y£ 3 L^4>N m - 

-Yg L^4>N m - Y£ L^4>N 3R + h.c, (11) 

with = 2020*. Using the triplet scalar, we can also have a Majorana mass term for the 
left-handed neutrinos: 

Anass = -Y? (W 7 (^2 A) - Y£ T Jj^f (ia 2 A) L tL + h.c. (12) 

In total, considering Eqs. (10), (11), and (12), the neutrino mass term can be written as: 

C maaa = ~^M^ + h.c, (13) 
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with * = {(u eL f, (u^f, (v t l) C ,N 1r ,N 2R ,N 3R ) T and 
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where we have defined m D l 



VffXg 1 and = v^Y^ p . and v& are the vacuum 
expectation values (VEVs) of the scalar doublet <fi and the triplet A. 1 We can diagonalize 
the symmetric mass matrix Ai u with the help of an orthogonal matrix O. We can study 
three different regimes for the parameters of the mass matrix: 

• m D l <C rn£ <C (separation scenario), 

• m2 <C <C M]j (type II see-saw scenario), 
~ m D l <C (hybrid scenario). 
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For simplicity, throughout our discussion, we will assume the parameters m 
Mjl to be real. We will show in the following how soft-breaking terms could alter the 
explicit form of the eigenvalues associated with Aii and M. v . 

The mass term for the charged leptons will also experience restrictions by the J-"-symmetry 
and is given by: 
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which can be rewritten in matrix form: 



e R 

£mass = - (el, W> Tl) Ml I flR | + tl.C, 

with the mass matrix hA\ given by 
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The charged lepton mass matrix Aii is diagonalized by a bi-unitary transformation, 

f m e 



U\M l U R = Mf^ 





m T 



1 Note that the VEV of the latter is constrained by corrections to the p parameter to be at most about 
1 GeV. 
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where m e , m M , m T are the electron, the muon, and the tau masses, respectively. We can 



then assume the presence of terms that break the L e 
lAJ 7 ! = 2 Dirac mass terms: 



Ln — L T symmetry softly through 
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where we can define s^f = v<pYp P and a,(3 = e,fi,r. Considering, for simplicity, a sym- 
metric form of the charged lepton mass matrix (s 
have: 
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Under the assumption of small and m D , we obtain 
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If we identify s e ^ 
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m T , we get an AAiAA] 



similar to Eq. (9). In this case, we would have a charged lepton mixing matrix Ul like 
the one in Eq. (5). We will show in the following discussion that the neutrino mixing 
matrix is compatible with the bimaximal form, even in the presence of soft-breaking terms. 
Therefore, the PMNS matrix will lead to the values of the neutrino mixing angles reported 
in Eqs. (6), (7), and (8). We assume the problem of 9 12 to be solved in this way. The task 
to investigate the details of a corresponding flavon potential will be left to others, as it 
would lead away from our main point, which is the prediction for the neutrino sector. We 
want to stress that any model with L e — L^ — L T symmetry will lead to similar predictions 
concerning the neutrino mass spectra, no matter how the problem of getting the correct 
9 12 is solved. 



3.1 Case 1 (separation scenario): <C <C Mjj 



Under the assumption that <C <C M 1 ^, with i. 
obtain: 
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with A± = ± ^{mf) 2 + (mf ) 2 and A± = ± ^(M R 2 ) 2 + (M* 3 ) 2 . Since A_ and A_ are 
negative eigenvalues of the mass matrix M. v , we choose to define the neutrino mass basis 
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as x = (^l-R) ^3ij, —iv5R, v§r) t ■ The orthogonal matrix is then given by 
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where we have set 



m L1 , mf 



M R1 and Mi 3 



m L2 , M l R 2 

interaction basis \l/ is related to the neutrino mass basis x by 
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This form of the light neutrino mixing matrix U u is exactly bimaximal if 9 = | (cf. Eq. (4)), 
which will happen in the limit run = vn^ 2 . To avoid the presence of two massless eigen- 
states (7713 = and M3 = 0), we can assume that the L e — L M — L T flavour symmetry is 
softly broken by |AJ-"| = 2 Majorana mass terms: 
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In this case, the following mass matrix results: 
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where si 01 



v/±Y^ a and a = e, //, r. Considering for simplicity s 
-C -C Af]j, we obtain after diagonalization 
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with A' ± = s ± ^(m e L y + (mf ) 2 , A s = s, A' ± = S ± ^(Afi 2 ) 2 + (M* 3 ) 2 , and A, = S. In 
case of 3 < M]j 2 , M^, 3 we have A' ± ~ ± v / (Mj ? 2 ) 2 + (M R 3 ) 2 . In this way we could explain 
the presence of one keV sterile neutrino (S ~ keV) and two heavier sterile neutrinos 
{M R 2 , M R 3> S), as required by the working example reported in Ref. [21]. For the light 
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neutrinos, we have m x = s + b, m 2 = s — b, and m 3 = s, where b 
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Note that we need s < due to the condition \mi\ < {m?], which leads to physical light 
neutrino masses of \m\\ = b— \s\, jm^l = b-\-\s\, and \m^\ = \s\. The resulting light neutrino 



mass square differences, in case of |s| <C \m e ^ ,eT \, are given by Am| x ~ 4(— s)b > and 



Am 2 ! 



2b(—s)—b 2 < 0. In this model, we predict a neutrino mass spectrum with inverted 
hierarchy. To fit the experimental data, we need b = 0.0489 eV and s = —3.9 x 10~ 4 eV, 
which allows to predict the light neutrino mass spectrum explicitly: \mi\ = 0.0486 eV, 
\m 2 \ = 0.0494 eV, and |m 3 | = 0.0004 eV. 

No matter if soft breaking is present or not, the mixing from the neutrino part is in any 
case given by Eq. (25), which can be reformulated as 
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Bimaximal mixing is perfectly restored in the limit m^i = m^, which should approxi- 
mately be true. Considering then Eq. (5) for the charged lepton mixing matrix Ul, we will 
then obtain neutrino mixing angles compatible with the experimental values. 



3.2 Case 2 (type II seesaw scenario): <C <C Mjj? 

Under the hypothesis that ~ mo, the eigenvalues of the 6x6 neutrino mass ma- 
trix are S = {A + , A_, 0, A + , A_, 0}, with A± = ^ ^^^ ^ + 0[m%) and A± = 

±-^/ (M R 2 ) 2 + (M R 3 ) 2 + 0(m 2 D ). As in the previous case, in this scenario we need soft 
breaking terms to avoid the presence of two zero eigenvalues. In this case, we have to 
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diagonalize the following mass matrix: 
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Taking into account that the eigenvalues of the matrices Mr are much bigger than the 
entries of the matrices rriL and m^, we can block diagonalize the matrix M. u in the 
following form: 

^(2),biock __ ( m L — m D M R 3x3 
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ing the matrices rriL — moM^m^ and Mr separately, we obtain the eigenvalues £' 
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S ± v^Mr, and A s = 5. Also in this case we could explain the presence of one keV sterile 
neutrino (S ~ keV) and two heavier sterile neutrinos (Mr ^> S 1 ). The light neutrino mass 
matrix is given by 
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where we have neglected terms of the order O (^jp^j an d where we have assumed that 



b = rriL — j^- > 0. Similarly to the first case, we predict neutrino masses that are given 
by rrii = s + b, m>2 = s — b, and m 3 = s. In this case, we again have b = 0.0489 eV, 
s = -3.9 x lO" 4 eV, |mi| = 0.0486 eV, \m 2 \ = 0.0494 eV, and \m 3 \ = 0.0004 eV, just as in 
case 1. 
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The light neutrino mixing matrix is given by: 
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e = 0), we would end up with a mixing 
matrix that is exactly given in Eq. (4), but the soft breaking is able to alter the form of the 
matrix. It is hence worth to investigate if e can be small at all: We know from the condition 
that the mass spectra have to be compatible with all experiments and observations that 
S ~ C(keV) < M R , m D < 0(100 GeV), and that (m L -jfc) = 0{m u ) > s. Furthermore, 
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Here, the first factor S/Mr is always small. The second one has a denominator mi — ^p., 

2 

which has to be of the order of the neutrino mass. In the limit itll ^> j^, this factor will 



be roughly equal to 



171 D I Mr 



and hence tiny. Even in the limit thl <C jf- it can at most 



m L > — - m r 

be of order one, which is not enough to compensate the smallness of the first factor. Only 

2 

in the extremely fine-tuned case, m v <C jjf-, it may be that this factor gets sizable, 
but this case is not to be expected. Considering the expression reported in Eq. (5) for 
the charged lepton mixing matrix Ul, we find neutrino mixing angles compatible with the 
experimental values. Accordingly, this model provides a natural benchmark scenario for 
keV neutrino Dark Matter that even works without additional requirements such as the 
one needed for case 1. 
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Using the same approximations as for the scenario in case 2, we obtain the eigenvalues 



S = {A' + ,A'_, X 8 , A' + , A'_, A s }, with A' ± = s± y/2m L ( 

S ± V2M R , and A s = S. Similar to before we have b = 0.0489 eV, s = -3.9 x 10~ 4 eV, 

|mi| = 0.0486 eV, \m 2 \ = 0.0494 eV, and \m 3 \ = 0.0004 eV, where b = \f2m L (l - jfe). 

Note that (due to < M#) b > is the only possibility in this case. 
The mixing matrix for light neutrinos is given by 
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In this case, e will always be tiny and could even be taken to be zero. Therefore, also this 
model will lead to neutrino mixing angles compatible with the experimental values. 



4 Conclusions 

We have studied how a L e —L^—L T flavour symmetry can be used to simultaneously explain 
observed masses and mixings of the light (active) neutrinos and how it simultaneously can 
lead to a heavy (sterile) neutrino sector, where one state is light. Such a flavour symmetry 
would then be responsible for a scenario with keV sterile neutrinos, which are a prime 
candidate for Warm Dark Matter. The L e — — L T flavour symmetry chosen actually 
predicts two neutrinos with nearly the same mass and one which is massless, in the sterile 
as well as in the active sector. Soft symmetry breaking, however, creates a non-zero mass 
for these particles, which will be much smaller than the masses of the respective degenerate 
pairs. Furthermore, the symmetry predicts zero #13 and maximal 6*23 and 6 , 12 , where the 
last prediction has to be corrected by the mixing coming from the charged lepton sector 
in order to lead to the leptonic mixing we expect to be present. 
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